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ABSTRACT 


Thermal images of a ship target were recorded with an AGA 
Thermovision 780 thermal imaging system on 6-13 May 1987. 
These images were used to produce temperature distributions 
ef the ship. The temperatures predicted with the AGA 
computer-processing algorithm were Gempared to actual 
temperatures measured with thermistors at eight locations on 
the ship's superstructure. An empirical modification based 
on the identity of the ambient and local atmospheric 
conditions was developed for the atmospheric compensation 
aigoritnm. The predicted temperatures were found to agree 
with the actual temperatures within 1.5 °C in 77 % of the 
measurements. The Contrast Transfer HMNeGEMem  {CTE), 
Modulation Transfer Function (MTF), and Minimum Resolvable 
Temperature Difference (MRTD) were determined for the 
Thermovision using flat black painted aluminum bar targets. 
The resultant curves showed the expected form with some 


experimental scatter at higher spatial frequencies. 
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ifn Le ODUCTION 


Thermal imaging is the process of collecting thermal 
radiation emitted by an object and producing a visible image 
of that object. Although the concept is simple, 
implementation of a thermal imaging device quickly becomes a 
complicated venture. First, the designer is confronted with 
a relatively broad portion of the electromagnetic spectrum, 
the infrared region; therefore, he must consider target 
radiation characteristics in the design process. Second, the 
target must be recognized against a background of several 
natural radiation sources. These sources include celestial 
emissions, the sun, the sky, and terrestrial materials. 
Finally, the atmosphere attenuates target radiation by 
scattering and absorption mechanisms. Most current thermal 
imagers utilize quantum detectors uo collect infrared 
radiation from a scene and convert this information to analog 
electrical signals which are processed and displayed ona 
video monitor. The ability of a detector to recognize a 
target depends on image contrast, i.e. the difference between 
target emittance and the background emittance. Because 
thermal radiation is a function of an emitter's temperature, 
the effective temperature difference between target and 
background provides a good measure of image contrast provided 


the emissivity differences are not great. 


The apparent temperature perceived by an imaging device 
can vary considerably from the actual temperature. Thi. 
irregularity is primarily due to background and atmospheric 
effects. If a designer is able to predict the impact of 
these effects accurately, then the imaging device will be 
better able to match perceived target temperatures with 
actual target temperatures. 

This thesis explored the differences between the apparent 
temperatures sensed by a thermal imaging system and the 
actual temperatures of a target. The imaging system used was 
an AGA Thermovision 71306 manufactured by the AGEMA 
COnporacion: The target was the R/V Point Sur, an 
oceanographic ship operated by the National Science 
Foundation for the Naval Postgraduate School and Moss Landing 
Marine Laboratory. The background was Monterey Bay in two 
vicinities: the Moss Landing Marine Laboratory and the 
Hopkins Marine Station. Data were taken during May 1987. 
Also, the Thermovision's Modulation Transfer Function (MTF) 
and Minimum Resolvable Temperature Difference (MRTD) curves 
were determined using square aluminum plates painted with 
Dawe sO “Tt labyblackeoatne. The use of aluminum permitted 
adequate heat transter “aeuecs a plate, so that the 
temperature was constant at all portions of the plate. The 
emissivities of the paint and aluminum were substantially 
different which caused a radiation contrast to exist between 


the bars of paint and bars of exposed aluminum. 


Chapter II details the theory behind thermal radiation, 
atmospheric ~ propagation of such radiation, and thermal 
imaging system performance measures. Chapter III describes 
the technical characteristics and operational capabilities of 
the Thermovision 780, as well as the temperature measurement 
techniques of the system. Chapter IV contains the results of 
the experiment, while the analysis of these results is 
presented in Chapter V. The conclusions and recommendations 


of this thesis follow in Chapter VI. 


Il. PERTINENT Poeer 


A. THERMAL RADIATION LAWS 

Thermal radiation is comprised of photons that are 
emitted by molecules undergoing vibrational and rotational 
Rais leronce These molecules may be the original sources of 
the photons or they may be reradiating photons emitted from 
other sources. In macroscopic terms, objects can be self- 
emitting thermal radiation sources or reflectors of such 
radiation. Most objects in our environment exhibit both 
mechanisms. “{(Ref> 1:p.) ve) 

Several physical laws have been developed that describe 
the nature of thermal radiation. Foremost is Planck's 
radiation law which gives the spectral distribution of 
radiant emittance, W\ 2 Dhaswlancwc. 

Wy [A,T] = (2mhc/d?) 1 (Watt/m* - um} (2-1) 

[exp(eh7NkT) = 2) 
where k is Boltzmann's constant, c is the speed of light in a 
vacuum, and h is Planck's constant. Figure 2.1 {Renae 
Dp... %23)] Sdepiecrs Spectral emittance for three object 
temperatures. When this expression is differentiated and set 


equal to zero, Wien's displacement law is obtained. This is: 
lee cones { um-K } (2-2) 
where ,, is the most abundant wavelength of radiation) fone 
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particular temperature. This thesis deals with the 8-14 
micrometer band Oreo 1 ac lom which corresponds coum 
temperature range of 207.0-362.2 K, temperatures that clearly 
exist in our environment. This is a good indication of the 
background problem that confronts thermal imaging systems. 
Integrating the expression for Planck's law across the 
entire spectrum of wavelengths yields the Stefan-Boltzmann 


law, given by: 
W(T) = (2m°k4)/(15c7h3) T4 = oT4 8 (Watt/m?} (=n) 


where W(T) is radiant emittance and o is the Stefan-Boltzmann 
Samscant . Radiation sources are commonly described in terms 
of radiance, L(T), which is radiant emittance through a unit 
solid angle. Many sources effectively radiate into a 
hemisphere as perfectly rough planes. Radiance for such a 


Lambertian surface is: 
Mec) = WC Ly / 1 {watt/m?-steradian} (2-4) 


Most objects in our environment do not emit as ideal 
blackbodies, but emit only a fraction of the blackbody 
radiant power. This fraction is known as the “emissivity" 


and can range in value from zero (a perfect reflector) to one 


(a black body). Thus, the Stefan-Boltzmann law becomes: 
W(T) = e0T4 (2-5) 
where € is the emissivity of an object. For an object in 


i 


thermal equilibrium with its surroundings, the objeeume 
absorbed power must equal the emitted power. This leads to 


Kirener.s Ss lanewhwenerse 
Sy CLA) (Zoe 


where a(A) is the object's spectral absorptivity. 

Thermal radiation impinging on an object can undergo one 
of three mechanisms. First, the radiation can be absi@meed 
and subsequently reradiated as described by Kirchoff's law. 
Second, the radiation can be reflected off the objec 
Third, the radiation can be transmitted through the Gbgeem 
Hence, the boundary conditions at an object's surface are 


described by the total power law: 
a(XA) + p(A) + T(A) = 1 (2=7) 


where p(A) is the spectral reflectivity and T(X)} is the 


spectral transmissivity. 


B. NATURAL RADIATION SOURCES 

Celestial background radiation includes infrared emission 
from solid material within our galaxy, as well as from 
extragalactic sources. Irradiance levels from celestial 
sources have been found to be around 107+6 W/cm?-um for the 3 
to 26 micrometer region [{Ref. 2:p2 3-238 Such levels are 
several orders of magnitude lower than those of other 
sources; therefore, celestial sources will not be considered 
in this thesis. 


eZ 


The sun can be treated approximately as a blackbody at 
5900 K, when the earth's atmospheric effects are ignored. 
Picunem2 2a RCE. 21p..  s=34 | depicts the spectral irradiance 
of the sun outside and within the earth's atmosphere. At 
5900 K the most abundant wavelengths radiating from the sun 
are in the visible region; however, a small percentage of the 
sun's spectral emittance is in the infrared region i.e., 
0.0986 percent in the 8-14 micrometer band, which equates to 
iowomeW7ms so NfPRef. 2:p.3-36]. Contributions of the sun to 
background radiation must be considered during daylight hours 
and particularly in background scenes that approach the sun's 
direction. 

The sky contributes to the background radiation through 
the scattering On solar radiation and emission from 
atmospheric constituents [{Ref. 2:p. 3-71]. Figure 2.3 [Ref. 
2:p. 3-71] depicts these contributions as well as solar 
radiation and radiation from a cloud. The presence of clouds 
affects both the scattering of near-infrared solar radiation 
and the sky's thermal radiation. Forward scattering in 
clouds is prevalent in the near-infrared region; however, 
multiple scattering dominates the forward scattering effects 
in a heavily overcast sky. The spectral emittance due to the 
sky emissions approaches that of a black body at atmospheric 
temperature in bands of high atmospheric absorption. 

Optically thick clouds are also well represented by a 


black body curve based on a cloud's average temperature; 


ILS 


however, the "wings" of this curve correspond to regions of 
high atmospheric absorption. Accordingly, the black body 
curve in these "wings" is based on the atmospheric 
temperature rather than the cloud's temperature. Figure 2.4 
(Ref. 2:p. 3-73] illustrates this effect of a cumulus cloud 
radiating as a blackbody at one temperature and the sky 
radiating as a black body at different temperature. 

While the radiation of numerous’ terrestrial materials 
could be examined, only that of marine backgrounds will be 
considered in this thesis. The marine background is affected 
by four factors: the optical properties of water; the surface 
geometry and wave slope distribution; surface temperature 
distribution; and bottom material properties (Ref. 2:p. 3- 
Osa 

Sea water has very low transmittance and reflectance 
values in the infrared region. The exception is that high 
reflectances do exist at large angles of incidence. When the 
sea surface is roughened by wind, the reflectance is 
Significantly reduced at angles of incidence near the 
AGI Zen. Figure 2.5 [Ref. 2:p. 3-106] depicts reflective, 
and emissivity of smooth water versus angle of incidence. 

As with any thermal emitter, the sea surface temperature 
determines the radiant emittance of the sea. The temperature 
gradient is quite pronounced in the upper 1.0 mm of the 
surface due to the cooling effects of evaporation. Figure 


2.6 [Ref. 2:p. 3-109] depicts this gradient. The gradivemaam 
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much less at lower depths but can be influenced by the 
convective activity of the water. Surface temperatures can 
be altered by surface contamination which restricts the flow 
of heat from lower depths; thus the apparent surface 
temperature will be lower than the bulk water temperature 
fret. 2:p. 3-109}. 

Because sea water has such low transmittance values for 
the infrared region, the properties of bottom materials have 


negligible effects on the background radiation. 


C. ATMOSPHERIC PROPAGATION OF THERMAL RADIATION 

The transmittance of the earth's atmosphere is less than 
one, since the atmosphere does not behave as a perfect 
dielectric. This transmittance is given by the Lambert-Beer 


law: 
Ta(A) = exp[-uU(A)R] C20) 


meme  tn(X) is the extinction coefficient and R is the 
distance of propagation. The average transmittance for a 
particular bandwidth i, to dz is: 
Ta = i M2 exp[-p(r dir - 
ie eee | P[-H(A)R] (2-9) 
2-1 JX 
Extinction of thermal radiation is caused by absorption and 


scattering processes. Thus, the extinction coefficient is 


given by the expression: 


U(A) = K(A) + a(d) (2-10) 
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where k(XA) is the absorption coefficient and a(idA) is the 


Scattermignecer ficient. Both coefficients have components 
due to air molecules, as well as aerosol particles. 
Therefore: 

KON) = Ket) to kas) (2-11) 
and 

A(X) = An(A) + Aa(A) (2-12) 


where the subscripts m and a denote molecule and aerosol 
respectively. The values of these coefficients have been 
empirically determined and depend on the density and 
composition of both the molecules and aerosols. 

Molecular absorption is due primarily to water, carbon 
dioxide, ozone, nitrous oxide, carbon monoxide, and methane 
(Ret s.-29 OD = “5=1012= 255-105). These absorbers limit thermal 
radiation in the atmosphere to two windows: the 3.5-5.0 and 
8-14 micrometer bands. 

Aerosols have been categorized in four standard types. 
Maritime aerosols are made up primarily of salt particles 
which act as condensation centers fOr water. The 
concentration of these particles is largely dependent on wind 
speed, while the size distribution is influenced by relative 
humidity: Continental aerosols consist of silicon, iron, 
sulphates, and organic material. One third of these aerosols 


can act as condensation centers. Urban aerosols are composed 
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of combustion and industrial PeoOgucies. PL eaeosphier ic 
aerosols contain Sulphate with an occasional addition of 
OUCenbe rust. [RET we Seppe 17-14 | 

Scattering by molecules in the 8-14 micrometer band is 
relatively insignificant. Wavelengths in this region are 
much larger than particle size, so that Rayleigh scattering 
would be valid. However, the X\~* dependence of Rayleigh 
scattering reduces the effect of molecular scattering. Mie 
scattering theory Pestana eo BeOciameleleS SOQ dll Sizes ; 
therefore, it is valid for aerosols. Because Mie scattering 
is somewhat independent of wavelength in the large particle 
limit, aerosol scattering is quite significant in the 8-14 
micrometer band. Mie scattering is most effective around 
r/X ~ 1, while particle size distributions extend out to tens 
of micrometers. 

Due to the complex nature of atmospheric propagation of 
thermal i 2(olsl she Lolo several computer models have been 
developed to predict atmospheric transmittance. LOWTRAN 6 is 
the model which will be used in this thesis. LOWTRAN is a 
FORTRAN based code that calculates atmospheric transmittance 
and radiance for a specified range of wavelengths along a 
designated path length. Seven atmospheric models are 
available including the U.S. Standard Atmosphere and a model 


based on user specified meteorological data. 
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D. PERFORMANCE MEASURES OF THERMAL IMAGING SYSTEMS 

A fundamental parameter restricting performance of any 
imaging system is contrast. In the case of thermal imagers, 
radiation contrast is used to assess how well a target can be 
seen against its background. The expression for radiation 
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where Wp is the target radiant emittance and Wp is background 
radiant emittance. Figure 2.7 [{Ref. l:p. 29] shows radar 
contrast curves for four background temperatures. 

The radiation contrast of bar targets can be used to 
construct a Contrast Transfer Funceron —(e05). As a bar 
target represents a square wave input to an imaging system, 
the CTF describes the system's square wave amplitude response 
at the spatial frequency of that target [Ref. 4:p. 114]. The 
CTF of a system is determined by plotting the radiation 
contrast for a range of spatial frequencies. The resultant 
function has an initial value of 1.0 (100 percent contrasm 
at zero spatial frequency and drops to a final value of 0.0 
(mo contrast) at the system's cutoff frequency, i.e. the 
maximum spatial frequency above which a system can no longer 
resolve a test target. 

A more useful measure of a system's resolution capability 
is the Modulation Transfer Function (MTF). "The MTF is the 


Sine-wave spatial frequency response" [Ref. 4:p. 114]. The 


ir 


MTF can be determined from the CTF using the following 


relationship [Ref. 4:p.117]: 
M(£) ~ (m/4)[C(£) + C(3£)/3 - C(5£)/5 + C(7£)/7] (2-14) 


where M(f) is the value of the MTF at ae single spatial 
frequency, f. C(f) is the value of the CTF at that 
frequency. BPreame ceo = |RefL. Typ. I9)]™econtains an MTF curve 
for an example imaging system. 

While the MTF is a valuable measure, it describes only a 
small part of an imaging system's ability to perform intended 
tasks. Ideally, a performance measure must be based on 
fundamental system parameters and "must relate to the 
performance of the system as it is intended." The most 
widely used performance measure for infrared detection 
systems is Noise Equivalent Temperature Difference (NETD). 
NETD iS a measure of ae system's ability to detect small 
Signals in noise. NETD is the target-to-background 
temperature difference that corresponds to a system's peak 
Seronal tO rms noise ratio equal to one. ([Ref. l:p. 166] 

NETD is derived from an expression for spectral irradiant 


power received by the detector. Such an expression would be: 
Px(A,T) = [W\(A,T)/] OAT, (Watt/um} (2-15) 


where W),(A,T) is the target's spectral emmittance, Q is the 
Selawancglessubcenaed by the system's optical aperture, Am is 


the target area from which radiation is received by the 
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system's optical aperature, and T, is the “system opuream 
Cransmission. Since the solid angilie@rs Ao/R and the target 


area iS aBR@, equation 2-15 can be written as: 
P, = (W)/T) AoaBT, (2-16) 


where A, is the area of the optical aperature and a and 8 are 
the system's horizontal and vertical subtense angles. 

The target's differential change in irradiant power with 
respect to temperature is of primary interest. By 
differentiating equation 2-16 with respect to temperature, 


such an expression is found: 
SP) /0T = [(GbAgi on) onw om (2-17) 


The system's differential change in signal voltage is 
found by multiplying equation 2-17 by the system's 
responsivity, R(A), a parameter that gives the ratio of 


Signal voltage output to incident power: 

aVe/dT = R(X) [(ABAgT,)/m] OW) /dT (2-18) 
Responsivity is given by: 

R(X) = [Vp_D*(dX)] / (abB,)}/2 (2=1o8 


where V, is the detector rms noise voltage produced by By, 
the noise bandwidth of a test reference filter, a and b are 


detector dimensions, and D*(X) is the system's specific 
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Gereewn ly mince LeI ng sequation 2-19 into 2-18 yields: 
aV5/8T = D*(X) (GBAgT,G) / [1(abB,)1/2] awy/atT [2-203 


Assuming small target-to-background temperatures, the 
temperature differential oi target emittance can be 
approximated by the derivative of Planck's law at the 


Packoround temperature Tp: 
aW,/8T ~ (Co/ATR*) Wy (TR) (2-21) 


mlere C5 = 1.4388°% 104 um-K. Specific detectivity can also 


be written as: 


Dee mach (AN Boje SON 


= (2-22) 


= > 
0 BOr “Nee Ap 


where Se) is the wavelength corresponding to the peak value of 
system detectivity. Inserting equations 2-21 and 2-22 into 
equation 2-20 and integrating across the effective waveband 
Xz to se yields: 

awe om =) GbAct ov, Bhp) C5 AP Wy) (Tp) dA (2-23) 

T (abB, O72 - EO eg 

Assuming a small signal approximation and rearranging 
equation 2-23 produces an expression for the signal-to-noise 
mat LO: 


SNR = AVs/Vq = AT GBAGTG, , D¥ (Ap) Co AP Wy (Tp) dA (2-24) 
n(abB, j1 1/2 Ap “Tae AT 


JM 


Based on the definition of NETD, the SNR is Set equaieees 
one and AT becomes NETD. Thus, NETD in final form forme 
Scanning imager is: 


NETD = m(abBn)!/2 AY Tp? 


_T i Wy (Tg) Gd>7-2 (K} (2-25) 
apt, Bis Nn) C2 e 


Pa 

Although NETD is a satisfactory measure of a system's 
target detection capability, it is of limited utility for 
thermal imagers because this measure does not account for 
image quality. A more appropriate performance measure for 
these devices is Minimum Resolvable Temperature Difference 
(MRTD) which incorporates the system MTF and provides a 
measure of how well an imaging system can resolve a target. 
MRTD is defined as the blackbody target-to-background 
temperature difference in a standard test pattern at which an 
observer viewing the system display can resolve a target [Ref 
ep O wo Ue 1 oe 

An analytical expression for MRTD is based on a system's 
Signal-to-noise ratio in the image of one bar. This ratio 


ees 
SNR; = 4/1 M(f£) [AT/(NETD] (pl/2)71 (2-269 


where M(f) is the system modulation of the image (i.e. MTF), 
the AT/NETD factor is the electronic SNR measured at a NETD 
reference filter with a target-to-background temperature 
difference AT, and pif2 is the ratio of actual Sysuem 
bandwidth to reference bandwidth. 
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As the observor views the target, this perceived SNR is 
modified by ~ four factors. First, the eye operates on a mean 
signal rather than peak signal; therefore, the SNR is reduced 
Dye Lactor Of 2/71. Se conc cho ws ONRe ls impEevedsby ‘a 
temporal integration factor of (TaF) 1/2 Whee ers 1S the 
effective eye integration time and F is the system frame 
rate. Third, the SNR is improved by the eye's spatial 
integration of the bar height which is seven times greater 
than the bar width. This factor is (7/2 £p)1/2. Finally, 
the bandwidth ratio is improved by the eye's matched filter 
Bec lona ei Ret. Ip. 186] 

A good approximate value of the perceived SNR has been 
found to be 4.5 for a 90% probability of detection {Ref. l:p. 
ioe |) - Using this value and the observer's perception 


factors, equation 2-26 is rearranged so that: 


MRTD(f) = AT = 3 NETD p,1/2 eg Deus 
M(f£) (TAF) 1/2 


where o,l/2 is the improved bandwidth ratio due to the eye's 
matched filter action. Because the displayed noise in many 
systems is white [{Ref. l:p. 189}, p,t/2 can be given by the 


simple expression: 
p,t/2 = (4af/n)1/2 Ts) 


Inserting equation 2-28 into 2-27, the final expression 


for MRTD becomes: 
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MRTD(f) = 6 NETD f(ap)1/2 (K} (2-25 
TGs Cuma 2 
Figure 2.8 contains an MRTD curve for an example imaging 
system. 

The NETD listed for the AGA Thermovision 780° is Ovum 
at 22 mee (Refs Sie The MRTD is not provided by the 
manufacturer but was determined during the experimental 
portion of this thesis. These results are found in Chapter 
Ev The MRTD deduced is then applied to the discussion of 


the measured temperature distributions in Chapter V. 
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III. DESCRIPTION OF AGA THERMOVISION 780 


A. GENERAL 

The AGA Thermovision 780 thermal imaging system 
incorporates a dual scanner which senses thermal radiation in 
two spectral bands and produces electronic video signals. 
These signals are amplified and used to display images on a 
black and white monitor. The scanner and monitor are shown 
mmeeigure 3sde [Ref .6)). The system is augmented with a 
microcomputer that displays digitized pictures of the 
system's video on a color monitor and records these pictures 
emeo a floppy disk. The computer uses a software program 
that color codes each picture according to the scene's 
temperature ehis Cebu tien . A schematic of the total 


Thermovision system is displayed in Figure 3.2 [Ref. 7]. 


B. DUAL SCANNER 
The dual scanner is actually comprised of two scanning 

systems. The shortwave system, 3 to 5.6 micrometer, has a 
Single Indium Antimonide photovoltaic detector and silicon 
lens. The longwave system, 8 to 14 micrometer, has a single 
Mercury Cadmium Telluride detector and germanium lens. Each 
system has av7° xX 7° field of view with a standard lens. A 
3.5° X 3.5° lens is also available for the long wave system. 
Each system uses vertical and horizontal scanning prisms to 


produce a 4 to 1 interlaced raster scan. Each of the four 


Z5 


fields has 100 scanning lines, although only 70 lines ame 
used to produce images. The interlaced frame in this 
configuration consists of 280 scanning lines. The detectors 
are cooled to 77 K with Dewar flasks containing liquid 


nlerogen: 


C. BLACK AND WHITE MONITOR CHASSIS 

The dual scanner is connected to two black and white 
monitors, so that scenes in the two spectral bands can be 
viewed simultaneously. The monitor chassis contains controls 
for adjusting the brightness and contrast of the display. 
The chassis also contains controls for adjusting the thermal 
level and thermal range of the system. These two adjustments 
are measured in Thermal Units which are arbitrary units of 
measure proportional to the intensity of the system's 
received thermal radiation. The thermal level control 
adjusts the DC level of the AC video signal, while the 
thermal range control limits the dynamic range of this 
Signai, 

These controls can be used in conjunction with a 
GallibEeariton DIOL to assess manually the temperature 
distribution of a Scene. The relationship between thermal 
level and thermal range settings versus temperature is shown 
in the calibration curve enclosed as Figure 3.3 [Ref. 6}. A 
change in thermal level corresponds to a nonlinear change in 
scene temperatures the Thermovision can sense. An adjustment 
of the thermal range determines the measurable temperature 
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range which iS approximately centered about the "median" 


temperature established by the thermal level setting. 


D. IF 800 MICROCOMPUTER 

A BMC IF 200 microcomputer processes the infrared data 
for the AGA Thermovision 780 and provides an automated means 
of assessing temperature distributions. On command, the 
computer's DISCO 3.0 program (AGEMA Corporation proprietary 
software) creates an eight color picture based on thermal 
values provided by the system and parameters inserted by the 
user. The color scheme depicts the temperature distribution 
ae 2 scene with each color representing a particular 
temperature range. The program provides temperature data to 
the user in several formats. The feature utilized in this 
thesis is the ability to assess the temperature associated 
with each pixel in the display screen comprised of an array 
of 128 X 128 pixels. 

The computer in this temperature evaluation mode displays 
a crosshair which is moved about the screen by means of four 
cursor control keys. The user positions the crosshair over a 
particular pixel and the computer displays the pixel's 


temperature in the lower right corner of the screen. 


E. TEMPERATURE MEASUREMENT TECHNIQUE 
The thermal measurement technique utilized by the 


Thermovision 780 is based on the relation: 


eee Pa tbat l-2,)Pe + (1-T,)P, ~ {Watt} ee ) 
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where P; is the total radiant power received by the system, 
Powis Gie radiant power from the object as a blackbody, P., is 
the radiant power from the object's surroundings as a 
blackbody, P, is the radiant power from the atmosphere waam 
blackbody, Tz is the atmospheric transmittance, and £4 ismieme 
object's emissivity. The emissivity of the surroundings and 
atmosphere are assumed to be one. 

The first term on the right hand side of equation 3-1 is 
the received radiant power emitted by the object. The second 
term is the received radiant power emitted by the object's 
Surroundings and reflected by the object. The third term is 
received radiant power emitted by the atmosphere. 

Because the system's thermal value measure is 
proportional to received radiant power, equation 3-1 can be 


WELEEeCO Nas . 
Ii; = Tafolo + Tall-Eg)Ig + (1-Ta)Ig (Thermal Units) mee 


where the I terms represent the thermal values of 
corresponding radiation sources. The received thermal value 


consists of two terms such that: 
1 ee ee (3-3) 


where L is the thermal level setting on the monitor chassis 
and iis a fractional portion of the thermal range. 
Substituting equation 3-3 into 3-2 and manipulating terms 


yields the following expression for the object*s thermal 
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Cy tee aes 7e4-1)l. - 1/€5(1/T,-1)Ia 138) 


The thermal values in the above expression are dependent on 
corresponding temperatures. The thermal value-temperature 
relationship is calibrated under laboratory conditions using 


the equation: 
I =A / [C exp(B/T)-1] (3-5) 


Meese A,b, and © are Calibration constants. Lice esteeeant 
Galidracion Curves are Similar to that shown in Figure 3.3. 
Equation 3-5 is used by the system to convert atmospheric and 
ambient temperatures to corresponding thermal values and to 
convert I to an object temperature. 

LiemenooDneslC Transmission factor is approximated by 


Ene system using: 
Ta = exp[-a(/d-1) ] (3-6) 


where a@ is an atmospheric attenuation constant and d is the 
Sesecance CEO object. The value of a specified by the 
manufacturer is 0.008 for the long waveband. A more accurate 
transmission factor can be computed using the LOWTRAN 
propagation code which is based on a form of Beer's law as 


given in equation 2-8: 


Ta(A) = exp[-H(A)R] (Se) 


a 


where wu(d) is the atmospheric extinction factor that 
represents scattering and absorption effects as discussed in 


Chapter II. Values based on this equation will be used in 


this thesis. 


30 


iVeeeacCOU ERED DATA 


A. RECORDING OF SHIP IMAGES 

The AGA Thermovision 780 was used to record thermal 
images in the 8-14 micrometer band of the R/V Point Sur at 
the Moss Landing location on 6 May and 13 May 1987 and at the 
Hopkins Marine Station on 9 and 11 May 1987. The map in 
Figure 4.1 depicts these locations on Monterey Bay and the 
approximate direction in which the Thermovision was aimed. 
Images were recorded with the ship presenting various 
aspects; however, only images of the ship at a 90° aspect 
were evaluated iol ee BS thesis. The Point Sur was 
instrumented with eight thermistors Ie measure the 
temperatures of various Parcs of the ship. These 
temperatures would be used as a basis of comparison for 
temperatures sensed by the Thermovision. Weather balloons 
were launched prior to each recording period. The radiosonde 
data recorded by the balloons were used as input to the 
LOWTRAN 6 program to determine atmospheric transmittance 
values. These data are contained in Table l. 

Images were recorded with the Thermovision's thermal 
range setting on two, five, and 10 Thermal Units; however, 
only those images recorded at the setting of five were used 
to evaluate temperature distributions of the ship. A thermal 
range of two Thermal Units did not correspond to a large 
enough temperature range in order to measure the entire ship. 
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A thermal range of 10 Thermal Units produced images with 
poorer temperature resolution than those images recorded at a 


setting of five Thermal Units. 


B. EVALUATION OF SHIP IMAGES 

Before the recorded images could be evaluated for 
temperature distributions, the Thermovision had to be 
initialized with the following values: the ambient air 
temperature around the ship; the atmospheric temperature; the 
emissivity of the ship; the range of the ship; and the 
atmospheric transmittance. The values of ambient temperature 
and the atmospheric temperature were assumed to be equal and 
were based on the air temperatures reported by ship 
personnel. These temperatures are contained in Table 2. The 
emissivity of the ship was assumed to be 0.95 as previously 
determined and reported [Ref. 9]. The range of the ship from 
the Thermovision's location was determined by measuring the 
length of the ship's image on the Thermovision's display 
screen and inserting Bhs value nbeliece: the following 


CrIGonometric relaueren: 
Ry =) Hoo cde) (S;/w)7+ = (wSo/tan 6) 1/8; {m} (4-1) 


where R is the ship's range, w is the display sceen width 
(l3cm), SQ is the ship length (41.2m), 6 1s Che SysStenaies 
of view (7°), and S; is the image length in cm. InSerGaamg 
values for the constants, equation 4-1 becomes: 

Ress See ze {m} (4-2) 


Sh¢% 


where the uncertainty in range was based on the differential 
method of error estimation. 

As mentioned previously, the transmittance values were 
calculated with the LOWTRAN 6 code. The atmospheric model 
was based on radiosonde data collected by the weather 
balloons. The Navy maritime aerosol model was used which 
Pacorpouated three aerosol components: a continental 
component; a "stationary" component produced by winds and 
waitecaps; and a "“fresh" component formed by current 
soma tions, (Refwme3:p. 17-35]. 

An air mass character of 2 was selected for the Moss 
Landing location, while an air mass character of 4 was 
selected for the Hopkins Marine Station. The air mass 
Character was a subjective rating of continental influence on 
maritime aerosols ona scale of one to 10. A rating of one 
would represent an open ocean, while ten would indicate 
strong continental influence. AS can be seen in the map 
enclosed as Figure 4.1, the Moss Landing location was quite 
exposed to the open ocean, while the Hopkins Marine Station 
was somewhat protected by the Monterey Peninsula. 

A horizontal path was specified with an altitude of 5 m. 
The range was based on the predicted range for a given image. 
Calculated transmittance values are contained in Table 2. 

The Thermovision's image evaluation program also had to 
be calibrated to a reference heat source prior to evaluation 


of the Point Sur images. The calibration was conducted on 
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13 October 1987 using a black body heat source adjusted to 
temperatures’ ranging from 6.0 °C to 30.0 °C. The black sgea 
consisted of a hollow aluminum cylinder covered with fiber- 
glass insulation. The cylinder was 28 cm long and 7.5 "cia 
diameter. A 2 cm entrance hole was located in the freomemer 
the cylinder. The rear surface was canted, so that incident 
light normal to the entrance hole would not be reflected out 
of the cylinder. A thermocouple was attached to the rear of 
the cylinder for the purpose of measuring the black body 
temperature. The interior of the cylinder was painted flat 
black. 

The calibration was conducted in a laboratory with a room 
temperature of 18.8 °C. In order to cool the black beavers 
below room temperature, liquid nitrogen was poured into the 
cylinder until the black body temperature was 2 °C. As the 
black body warmed to room temperature, the Thermovision was 
used to record images of the black body at 1 °C increments 
Starting at 6 °C. In order to warm the black body above roam 
temperature , the cylinder was wrapped with a nichrome wire. 
The wire ends were connected to a variac, so that the current 
generated would heat the cylinder. Again the Thermovision 
was used to record images of the black body in 1 °C 
increments up to 30 °C. 

Once the black body images were recorded for each 
temperature, the Thermovision's image evaluation program was 


used to predict the black body temperatures. The results of 
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PicseGaelorcationmeprocess are. contained in Table 3. These 
results were “wae to produce the calibration curve shown in 
Figure 4.2. This curve was subsequently used to convert 
temperatures computed with the Thermovision method to values 
that more closely approximated actual temperatures. 

The temperature distributions of the Point Sur were 
Geveloped using a 6 X 15 element array. ie Ese eee inaly 
configuration was based on the Point Sur's approximate height 
mee length ratio of 1/5. The vertical and horizontal 
dimensions were multiplied by factors of six and three 
respectively in order to facilitate the presentation of 
ILigteie gags 


icseeieuerons on 8 172" %X 11" notebook paper. 


mom 1gurdat1on one distribution element approximately 


represented a 2.75 m X 1.35 m portion of the ship's surface 
area at 90° aspect. Images used to produce these 


distributions were recorded during the following times: 


DATE TIME SLDE OF SHIP 
6 May eS 7 Star Dockd 

6 May 1901 Portside 

9 May E28 Portside 

9 May 1140-1141 Starboard 
11 May 0944-0948 Sherehaleleler cs! 
11 May 1004-1008 Portside 
13 May 52 O21 Starboard 
13 May 526 Portside 
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The temperature distributions for these periods are found 
in Figures 4.3 through 4.6. An anomoly was found during the 
initial review of these distributions. The inconsisteme, 
pertained to Figure 4.3 which portrayed the Point Sur's 
temperature distribution for 6 May. Some temperatures of the 
ship's upper superstructure in this figure were lower than 
the ambient air temperature. Physically, this’ sitW@eapien 
would not be possible. The Point Sur had no known heat sinks 
in this portion of the ship; therefore, ship temperatures 
could equal or exceed the ambient air temperature, but would 
not be less than the air temperature. The conclusion drawn 
from Figure 4.3 was that the Thermovision might have been 
underestimating actual temperatures. This problem had to be 
resolved before any detailed analysis of the Point Sur's 


temperature distributions could be made. 


C. THERMISTOR DATA 

The thermistor resistance values were recorded by 
microcomputer throughout the time the Point Sur was at sea. 
The locations of the eight thermistors are indicated in 
Figure 4.7. After the cruise these values were converted to 
temperature values and averaged for the time periods 
corresponding to those periods for which the Thermovision was 
recording images. Actually, the thermistor averaging time 
was extended two minutes prior to and two minutes after the 
the image recording period. The averaged thermistor 
temperatures with standard deviations are found in Table 4. 
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These averaged temperatures were then compared with the 
temperatures- computed by the Thermovision for areas of the 
ship where the thermistors had been located. The temperature 
differences between the thermistor-measured temperatures and 
the Thermovision-computed temperatures are listed in Table 5. 
These differences were quite significant on 6 May and 13 May. 
The Thermovision underestimated the thermistor-measured 
temperatures by an average of 4.5 °C on 6 May and 4.7 °C on 
13 May. This situation reinforced the tentative conclusion 
ema t the Thermovision was underestimating aG ea 


temperatures. 


ee cle, MIF, AND MRTD CURVES 

Because a Minimum Resolvable Temperature Difference 
(MRTD) curve was not provided by the manufacturer, one was 
produced for the purposes of this thesis. As discussed in 
Chapter II, the MRTD curve provides a good measure of a 
system's resolution capability. Since the MRTD curve is 
based on a system Modulation Transfer Function (MTF), the 
Thermovision's MTF had to be determined. A system Contrast 
Transfer Ja Dokere sie) gl (CLE) was more readily attainable; 
therefore, the system CTF was determined and the MTF computed 
From this curve. 

PecOmerase  Iransiter Function (CTF) for the Thermovision 
was determined on 19 October 1987. The CTF was based on 
values of radiation contrast for targets of various spatial 
frequencies as described in Chapter II. The targets consisted 
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of 1/8 inch thick square aluminum plates with four barsmes 
flacgmpiack paint applied. The emissivity of the paint was 
estimated at 0.95 and the emissivity of the aluminum was 
estimated at 0.30 [Ref. 10]. Thus, a target appeared to the 
Thermovision as a standard seven bar target as depicted in 
Figure 4.8 [Ref. l:p. 76]. Bar widths ranged in size from 5 
mm to 120 mm. The targets were placed at distances ranging 
from 1mto4m. The following relation was used to generate 


28 various spatial frequencies: 
f= O57 220 {cycles/radian} (4-3) 


where f is spatial frequency, d is distance between target 
and Thermovision optics, and w is target bar width. 

With the laboratory room temperature at 18.8 °C the 
targets were heated to 25 °C. A thermocouple was bolted to 
the center of each target to measure target temperature. 
While the Thermovision was imaging a bar target, “an 
oscilloscope was used to meaSure the radiation contrast 
between bars. The oscilloscope was connected to a video 
output port of the anologue-to-digital convertor and set to 
display target signal voltage. The intensity profile in 
Figure 4.8 depicts how a bar target would appear on the 
oscilloscope display. The heights of the signal peak and 
minimum were measured and equation 2-8 used to calculate 


radlationecomerasite 


C= (Wr = Wp) / (Wr + Wp) (2-8) 
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in this case Wp, waS proportional to the emittance of the 
black bars “of paint and Wp was proportional to the emittance 
of the aluminum bars. (Note: The vertical distance on the 
oscilloscope display was proportional to signal voltage and 
thus target emittance. ) Since the target temperature and 
background temperature were equal, the radiation contrast 


could be written as: 


c= (EpoT 4 = €pOTp*) / (EpoTp4 oF EROT 4) 


= (Em = Ep)/ (Ep ot ER) (4-4) 


The values Gye jetetolaehealfeyal Cone nace sOre evar] OUsme Spatial 
frequencies are listed in Table 6. These values were 
normalized relative to 4 cycles/radian, the frequency at 
which the bar target consisted of one black bar and one 
aluminum bar. This target filled the system's entire field 
of view, hence this was the minimum spatial frequency that 
could be viewed by the Thermovision. The normalized 
radiation contrast values were then used to calculate the MTF 


values for each spatial frequency based on equation 2-14: 
lew 4) (Ci 2) + Ci 3f)/3 eC (Stiy/S + C(7f)77) (2-14) 


Orr and MTF values are found in Table 7. Finally, the MRTD 
was calculated for each spatial frequency using equation 


2-29: 


MRTD(f) = 6 NETD f(agp)1/2 (229) 
M(f) (nT,F)l/2 


ae, 


The NETD move elie Thermovision was listed by the 
manufacturer as 0.12 °C at a background temperature Of 2275s 
However, the MRTD calculated in this thesis was based on a 
background Cemperature of 18.8 °C; therefore, a correcraem 
factor was needed. This factor was readily computed using 


the expression for NETD found in equation 2-25: 


NETD = m(abB,)1/2 pie 2 ff Pw (T) dea (2-250 
aBAoTS D* ike \co : 


The factor dependent on background temperature was: 


CF(Tp) = Tp fe W, (Tp) daj7t (4255 
et 4 

The denominator of this factor was integrated from 4, = seme 

zO Ap = 10 wumwith a hand-held calculator using Simpson's 


rule. The corrected NETD was: 
NETD(18.8 °C) = CF(18.8 °C) NETD(22 °C) = 0.985 (0033 ee 
“Or. °C 


Thus, the NETD was essentially unchanged. 

The following Thermovision parameters were inserted into 
equation 2-29: 

NET D== 70212. <¢ 

oe c pa Ue iets (6: 

Te = 0.2 sec 


F = 6.25 frames/sec 
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Thus, the expression for MRTD became: 
MRTD = 4.00 X 1074 £/M(f) (°C} (4-6) 


MTF and MRTD values are found in Table 8. A plot of the 
system CTF, MTF and MRTD values are found in Figures 4.9, 


eo, and '4rrr. 
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V. ANALYSIS OF DATA 


A. GENERAL 
As discussed in Chapter IV, the data presented in 

Figure 4.3 and Table 5 showed that the Thermovision-sensed 
temperatures were substantially lower than the actual 
temperatures measured by the thermisters on the four days of 
data recording. The average Thermovision-thermistor 
temperature differences were -4.5 °C, -1.9 °C, -1.1 °C, and 
-4.7 °C respectively. An explanation was needed as to why 


the Thermovision was underestimating target temperatures. 


B. REVIEW OF THERMOVISION'S THERMAL MEASUREMENT PROCESS 
As detailed in Chapter III, the Thermovision computes 


object temperatures based on two equations: 

I; = Tafoly + Tal(l-Eg)Ig + (1-Ta)Ta (3-2) 
and 

L*= A [Cexp(B/T) See (3-5) 


The first term on the right hand side of equationmaee 
represents radiation emitted by the target. The middle term 
represents radiation from the target's ambient atmosphere 
that reflects off the target and propagates to the 
Thermovision. The Enume term represents atmospheric 


radiation received by the Thermovision. 
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Peel eelesic lemmviccamassimed to equal I, ,i.e. the 
atmosphere was treated as a continuum extending from the 
target's surroundings to the Thermovision. Thus equation 3-2 


becomes: 


Based on this assumption, equations 3-2 and 5-1 hada 
shortcoming. First, the middle term on the right hand side 
of these equations had an implicit value of atmospheric 
emittance equal to unity, i.e. the ambient atmosphere was 
treated as a black body. iMic-eairva tenn on thew@maight mand 
ace Or these equations had 1-T,j as the emittance of the 
atmosphere existing between the Thermovision and target. 
Since the ambient atmosphere and the atmosphere between the 
Thermovision and target were to be treated as one body, both 
portions would have to have the same emittance. This 
Eiirtance would be 1-T.. The original assumption that the 
ambient atmosphere emitted as a black body is invalid. Thus, 


equation 5-1 becomes: 


This equation was used to recalculate temperatures in the 
temperature distributions; however, the discrepancy was still 
not resolved. Micm==aoOLOdem==at this OOlntc was to fit 
empirically the thermistor temperature data to a modified 


feestone or eCqudtion 5-2. This required setting the first T, 
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term of the middle term on the right hand side of equation 
5-2 to one. Based on this change and a combination of like 


terms, equation 5-2 becomes: 
I; = TH8olo + (1=Ta)C2-fop se (See 


Theoretically, this modification could not be justi 
but it did result in good agreement between Thermovision- 
computed temperatures and thermister measured temperatures. 
The implication of equation 5-3 is that the atmosphere 
exhibited no attenuation of ambient radiation reflected from 


the target. 


C. REVISED TEMPERATURE DISTRIBUTIONS 

Ocak legal yzaliere| equations 5-3 and 35% temperature 
distributions of the Point Sur were recalculated with a 
Hewlett-Packard Model 97 programmable calculator. The 
revised distributions are enclosed as Figures 5.1 through 
se 

A number of observations about these distributions could 
be made. The warmest part of the ship was usually the smoke 
stack area which had temperatures clearly exceeding 20 °C in 
most cases. The one exception was the starboard side on 6 
May. Element B7 was the warmest part in this case (19.8 °C). 

Excluding the smoke stack, the warmest part of the ship 


varied from day to day. These areas are indicated below: 
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DATE LOCATION (element) SIDE TEMPERATURE 


6 May - B7 Starboard 19.8 °C 
9 May All Port 2054. °C 
11 May Be ,Bo Sigaebodrd ZO Bas oe 
13 May eZ Seanboara 20 gai we FC 


The higher temperature associated with row A (just above 
the water line) were most likely due to additional thermal 
radiation from the ship reflecting off the sea surface. The 
higher temperatures associated with row B corresponded to an 
enclosed portion of the Syglaljo) Apparently, heat was 
transferred from these enclosed spaces to the outer surface 
of the ship. The higher temperature at C2 was possibly due 
monsOolar glint. 

The coolest part of the ship also varied on a daily 


basis. These areas are noted below: 


DATE LOCATION (element) S1DE TEMPERATURE 
6 May Ce3 Starboard 1532 Ve 
9 May A2 Starboard Less Ae 
11 May AlS Pore rans °C 
13 May A2 Port 148.62 


Generally, these cooler portions were at extreme points 
Soecne ship's superstructure. Exposure to the open air or 
sea surface at these points probably resulted in greater heat 
dissipation than at other locations of the ship. 

A final note on the distributions is that the temperature 


difference between adjacent distribution elements rarely 
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exceeded the estimated deviation in temperature due to 
experimental error (1.5 °C) except around the smoke stagm 
however, the range of temperatures within the entire 
distribution, disregarding the smoke stack, clearly exceeded 
this error estimate in all cases. This range was 4.6 °C, 

6.6 °C, 6.5 °C, and §.9 °C respectively for each of the teu 


days. 


D. PRECISION ESTIMATES FOR THERMAL MEASUREMENTS 
The estimate of experimental error was based on the 
differential method using equations 5-1 and 3-5. The error 


in total thermal units due to errors cme, d&,, and dT waz 


dI; = d1/aT, dT, + a1/ae, deg + 3I/dT AT 


[Eolo + (£5-2)I5) GTa + [Talg > (Teo) ealeces 


+ [(1;?%B/AT*) exp(B/T) ] (5-4) 


Subsituting in estimated values dT, = 0-01, d&, = 00a 
dT = 0.2 °C and typical values for the I terms, the estimate 
of thermal value error was di =. 75eac.u, 

The subsequent error in temperature calculations due to 


thermal value error was: 
dT = aT/a1, aI, = [(T*/1,3) 7 (eee cy (5-5) 


Subsituting typical values for temperature and thermal units, 
the estimated error in temperature measurements was 


GU g lO a. 
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An additional error in temperature values was introduced 
aiieing. ee process of making the 6 X 15 element temperature 
distributions from images displayed on the Thermovision's 
computer CRT. Depending on the range of the Point Sur, two 
to four CRT pixels were averaged to establish a temperature 
for each element in the 6 X 15 arrays. Temperatures varied 
by as much as 1 °C between CRT pixels. The estimated error 
attributed to this averaging process was 0.5 °C. Thus the 


total estimated temperature error was: 


Total error = Measurement error + Averaging error 


aCe eee a Cea 5 °C 


E. REVIEW OF REVISED THERMOVISION VS. THERMISTOR 
TEMPERATURE DIFFERENCES 


The recalculated temperatures were also compared to the 
thermistor temperatures. These results are found in Table 9. 
Based on review of Table a; Thermovision-computed 
temperatures agreed with thermistor temperatures (within 
experimental error) in 20 out of 26 cases during the four 
days of measurements. Six additional cases could not be 
compared because temperatures exceeded the Thermovision's 
thermal range setting. 

In the six cases in which agreement did not exist, few 
conclusive trends about these anomolies could be established. 
These cases involved four different thermistors over the 


course of all four days. One notable observation was that in 
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five cases the Thermovision's predicted temperatures were 
much lower than the thermistor temperatures. Apparently the 
thermal radiation was partially shielded or severely 
attenuated in these cases. 

A point that must be considered 1S that each thermaamen 
was meaSuring an area of ship surface on the order of one 
Square centimeter while each element in the Thermovision's 
temperature distributions corresponded to an area on the 
order of one square meter. Thus, any localized heat sources 
on the ship may have been obscured by more dominant ambient 


Gonadttions. 


F. MTF AND MRTD CURVES 

One observation of the MTF/MRTD plots (Figures 4.10 and 
4.11) was the cutoff frequency, i.e. the spatial frequency at 
which the Thermovision could no longer resolve a target. 
This frequency was 400 cycles/radian; however, this value has 
little utility because it corresponds to an infinite MRTD. A 
more practical frequency was the critical frequency which was 
determined by extending a line through the approximately 
straight portion of the MTF curve to the intercept of the 
frequency axis. This critical frequency was 360 cycles/ 
radian which corresponded to an MRTD of 2.2 °C. Asiana 
illustrative example: the critical frequency would equate to 
a lateral distance of 1.4m for a target at 1000em. Spas 


the Thermovision could resolve target structures to within 
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1.4 m provided that the target to background temperature was 
ae leasta2. 26°C. 

The MRTD curve was also used to determine the MRTD of an 
element from the 6 X 15 element temperature distributions. 
Recalling that these elements represented 2.75 m width along 
the ship's length, the horizontal spatial frequency of an 
element is 182 cycles/radian at a target range of 1000 m. 
This frequency corresponds to an MRTD of 0.2 °C. An 
important observation is that this MRTD is much less than the 


uncertainty in temperature due to experimental error 


me 86°C). [ise coils pixel Sige did Nore ani eae ne 
experiment. However, pixel size could be a factor if the 
MRTD was comparable to the temperature uncertainty. The 


limiting element size in this case would be 1.6m based on an 
mwoeor i.5) °C. Thus, the Point Sur with a length of 41.2 m 
could be divided into a maximum of 25 horizontal elements for 
temperature distributions. 

A clarification at this point should be made with regard 
to vertical resolution. The system MRTD curve developed in 
this thesis was based on horizontal spatial frequencies; 
therefore, the analysis of vertical resolution was not 
possible in this thesis. 

The Thermovision's computer displayed the system's 
7° X 7° field of view with the array of 128 xX 128 pixels. 


This pixel size corresponded to a spatial frequency of 525 
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cycles/radian. Thus, the image video would be limited by the 
Thermovision's components preceding the computer display. 

A final observation of the MTF and MRTD plots was that 
these curves were very typical except for a portion of the 
MTF at the higher spatial frequencies. In this region the 
Slope of the curve for a typical MTF would be monotonically 
decreasing. This was not the case for the Thermovision's 
MTF; however, this was most likely due to data scatter, 
rather than the Thermovision having atypical performance 


characteristics. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. SUMMARY 

Thermal images of the Point Sur were recorded at two 
locations on Monterey Bay with the AGA Thermovision 780. The 
Thermovision's temperature evaluation program was calibrated 
under laboratory conditions with a black body source set at 
known temperatures. This program was used to establish 
temperature distributions of the Point Sur images. These 
temperature distributions consisted of 6 X 15 element arrays 
which portrayed the ship's superstructure. 

Temperatures from eight elements on these distributions 
were compared to temperatures measured by thermistors. These 
measurement devices had been installed during the recording 
of the ship's images and were located at eight points on the 
Point Sur that corresponded to the eight elements’ in the 
temperature distributions. A deduced correction was made to 
an equation used in the Thermovision's temperature evaluation 
program, so that Thermovision-computed temperatures agreed 
with thermistor-measured temperatures. 

System CTF, MTF, and MRTD curves were also produced. The 
MRTD curve was used to determine the resolution limit of the 
Thermovision and to compare this limit with the resolution 
ema ct oie the temperature distributions and the 


Thermovision's computer display. 


One 


B. CONCLUSIONS 

Initial comparison of target temperatures computed by the 
AGA Thermovision 780 and actual temperatures measured by 
thermisters revealed that the Thermovision underestimated 
actual temperatures on the four days on which data were 
recorded. Two modifications were made to an equation used by 
the Thermovision's computer to calculate target temperatures. 
The first modification was made based upon the assumption 
that the target's ambient atmosphere was’ essentially 
identical with the atmosphere between the target and 
Thermovision. This assumption had the implication that the 
ambient atmosphere's emissivity would be the same as the 
other portion of the atmosphere. The second modification was 
empirically determined and when used in conjunction with the 
ELESE = Meal rTedcion:, good agreement existed between the 
Thermovision and thermisters to within 1yS ae in the few 
cases in which agreement clearly did not exist, the 
Thermovision underestimated target temperatures. This was 
attributed to the masking of target radiation rather than any 
measurement error introduced by the Thermovision. 

The AGA Thermovision 780 MTF and MRTD curves produced 
from images of aluminum bar targets resembled the curves of a 
typical system depicted in Figure 2.8 except at the higher 
Spatial frequencies of the MTF curve. The irregularities in 
this portion of the curve were attributed to data scatter 


associated with the very small target signals in these 


Bye 


frequencies. The Themmovssion's GuLOTEMerrequency was 
determined to be 400 cycles/radian. The critical frequency 
was 360 cycles/radian at 2.2 °C MRTD. 

Another important spatial frequency was that which 
corresponded to the width of an element in the temperature 
distributions. This frequency was 182 cycles/radian at 


a2. C MRID. 


C. RECOMMENDATIONS 

While the Thermovision's temperature prediction 
capabilities proved satisfactory under various atmospheric 
conditions, further testing should be conducted to verify 
that changes to one of the system's temperature prediction 
equations are justified for other conditions. 

The use of 6 X 15 arrays to depict the Point Sur's 
temperature distribution was viable but could be improved. 
First, the number of elements in these arrays could be 
increased in order to reduce the error introduced when 
several CRT pixels are averaged to determine a temperature 
for one array element. The maximum number of horizontal 
pixels would be limited to 25. The Thermovision's MRTD for 
this pixel size would equal the experimental error in 
memperature (1.5 °C). Second, an automated means of 
extracting temperature information from the Thermovision's 
computer is needed. The manual process of recording 
temperatures from the CRT screen and calculating average 
cemperatures is tedious and time consuming. 


Sys: 


The determination of the system MTF and MRTD curve could 
be improved by using an amplifier to increase the gain of the 
video entering the oscilloscope. This action would increase 
the accuracy of measuring radiation contrast and reduce data 
scatter at the higher spatial frequencies where target 
Signals are quite small and difficult to measure. Also, the 
Thermovision should be evaluated using the MRTD test 
described in reference one and the results compared to a MRTD 
curve determined using the method presented in this thesis. 
Finally, MTF and MRTD curves should be developed for vertical 
Spatial frequencies, so that the Thermovision's resolution 
characteristics in the vertical direction can be assessed. 

Any further testing with the Thermovision should be done 
with the Thermal Level and Thermal Range settings adjusted so 
that the entire target distribution of a target can be 
determined. One limitation of this thesis has been the lack 
of temperature information around the Point Sur's smoke 
stack. This limitation was due to the Thermal Level being 
set too low, so that the highest temperatures on the ship's 


Superstructure could not be measured. 


54 


as APPENDIX A 


FIGURES 


BS o 
Tg = 280 K 





Figure 2.1 Planck's law for spectral radiant emittance 


= 


Solar Spectral Irradiance Outside Atmosphere 


Solar Spectral Irradiance at Sea Level (m = 1) 
Curve for Blackbody at 5900 K 
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Figure 2.2 Spectral distribution curves related to the sun 
(Shaded areas indicate absorption due to 
atmospheric constituants) 
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Figure 2.3 Contributions to background radiation ines 
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Figure 2.4 Spectral radiance of cumulus cloud 
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Figure 2.5 Reflectance and emissivity of water ( 2 to 15 um 
average) vs. angle of incidence 
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Figure 2.7 Radiation contrast for 8-14 um band 
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Figure 2.8 MTF and MRTD curves of thermal imaging system 
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Point Sur temperature distribution (°C) 
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Figure 4.8 Bar target pattern perceived by Thermovision 
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Figure 4.9 Thermovision system CTF curve 
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Figure 4.10 Thermovision system MTF curve 
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Figure 4.11 Thermovision system MRTD curve 
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APPENDIX B 


TABLES 


TABLE 1 


RADIOSONDE DATA 


DATE 6 May 9 May WES alt 13 May 
MEME OF LAUNCH EOS 0947 10Cc9 1008 
PaIMPERATURE (°C) io OL Ite iPerr on! 14.01 
mew POINT (°C) LWA ae al d leae oy It Ie Aaalkal ereoas 
meio iVE HUMIDITY (%) S220 SCBA fey Ih 99.01 Ses 
PRESSURE (mbar) GO Sai eee 1S od 00 9r 71 iC Oier Sil 
WIND DIRECTION (degrees) 238.31 5) Sispemel il. Jes Aba Spoil acca 
WIND SPEED (m/sec) lero Son eel sesh 
Sire tTUDE (m) 7a.2) 610) th sO heal g470 2 Sysco al 


a 
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TEMPERATURES AT SHIP AND TRANSMITTANCE VALUES 


DATE 6 May 
TIME reso-iea6 
AIR TEMP. (°C) 14.5 
SEA SURF. TEMP. (°C) 15.2 
SHIP RANGE (m) 800 
TRANSMITTANCE Onl 
DATE 11 May 
EIME 0342-0950 
AIR TEMP. (°C) an 
SEA SURE. LEMP (7 i. 4 
SHIP RANGE (m) 800 
TRANSMITTANCE Oe Ske 


6 May 
1859-1908 


14.5 


ll May 
OO ZS Rene 


Zen 


9 May 


[26-2 


Ls May 
158-1528 


14.0 


9 May 


1138-ies 


13 May 
1524232 


14.0 


TABLE 3 


THERMOVISION SUPPLEMENTAL CALIBRATION DATA 


BLACKBODY THERMOVISION MEASURED TEMPERATURE 
TEMPERATURE TEMPERATURE DIFFERENCE 
(°C) (°C) e) 
oeOms OL See Ol =O. 9 tae. 
fae0 65. - 0.9 
ain 18) ae - 0.8 
9.0 ors + 7.8% 
ligule, 9.1 - 0.9 
ed. 0 Or - 0.9 
Jeg) ee aes - 0.7 
3... 0 2S = 0255 
Aa Lee's) - 1.5* 
5200 13.6 - 1.4* 
16.0 ILS 5 (é - 0.4 
io. 0) be. 5 - 0.5 
jefe) G8) ee 0 G0 
19.0 19.1 + Orr 1 
ZO. 0 20 2.0 Oz.0 
2a 0 OES Oro 
E20 22k + QO.l 
23.0 2324 On 4 
24.0 24.4 + 0.4 
25). 0 25.4 + 0.4 
26.0 Zoo + 0.6 
27020 26 + 0.6 
Ze .0 239 + Q.7 
Z9)..0 29°23 + 0.3 
BC. 0 Silia2 + 1.2 


* Data point disregarded. 
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DATE 


TIME 


THERMISTOR 


THERMISTOR 


THERMEST OR 


THERMIESLOR 


THERMISTOR 


THERMISTOR 


THERMISTOR 


THERMISTOR 


DATE 


TIME 


THERMISTOR 


THERMES TOR 


THERMISTOR 


THERMISTOR 


THERMISTOR 


THERMISTOR 


THERMESTOR 


THERMESTOR 


TABLE 4 


THERMISTOR DATA Gao. 


6 May 
FS55— ess 
cee sje. 0) ILC 
ee SO lee 
21 PEOr 
aes) Sise8) - 1G, 
Le .7 LEO a7 
LT <2ZoEO SEZ 
24.05+0.34 


217 620.3 


PT May 

0942-0950 
L5.41+03a8> 
PS 196 i04 06 
16.0720 506 
Se oO 
Eoin 2 340. 19 
16% 8920..13 
S56 0250 co 


39). toe Oe 


6 May 
1859-1963 
PG. 3EEO bs 
16: 3 726r is 
ZO Saeeemel 2 
bi fer 4/22 0) 6 2G 
18 .64207, 29 
16. oae0r 1S 
24 323304.55 


PRs e210) 6 718 


JI May 

1002-1010 
bey 746 )oc(6) GIL 
1.5709 20 
15.9 2202016 
15. 6682070 
1439 510).22 
1629020 320 
34257207201 


38 2420 9 
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9 May 
LI2Z6—1139 
16 7550 —e 
D7. Peon 
17.4120 
lL]. Zac 
1S. 2220eec7) 
17 - 6 UaOrene 
34.6340.24 


332 2220R60 


t3 May, 

IS16-T2s26 
1373 50 eene 
19% 31 eOrece, 
18)65+0ieus 
207 SiatOrez0 
20.98 LOO 
20 Vat Oued: 
$97 50 Oeaas 


39 26 Oa,08) 7.6 


9 May 
1138-1lla@3 
15.5840 2a 
17 .00tGRas 
17.6702 
16.7220@sas 
18.6/+£0e2 
17 .22£02ae 
33.5820" 


28.16+0.40 


13 May 
1524-1526 
17.96+0Reme 
18.3340 
20 .23£0Re 
20 .10+0%ez 
27.42£03ge 
18.6120 23s 
41.36+10.59 


37. 34t0Ram 


TABLE 5 
TEMPERATURE DIFFERENCE (°C) 
(THERMOVISION TEMPERATURES-THERMISTOR TEMPERATURES ) 


DATE 6 May aerial ll May 13 May 
TRSRMISTOR 1 = 36 a ees + 1.8 ee eno 
MEERMISTOR Z - 3.4 a agate ee - 1.4 = Ove 
MaERMiSTOR 3 = 15) 4% ma. o = lee 2 8 
THERMISTOR 4 = gee — Jit - j1.4 - 3.9 
MASRMESTOR § =m Se - 2.0 == ene ae aye 
THERMISTOR 6 eS 2 ae. © = ele hei 
MAERMISTOR 7 ee. 2 * * * 

MHERMISTOR 8 a SS x x * 


* Not available. 
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TABLE 6 
RADIATION CONTRAST 
(Target Temperature = 25.0 °C, Room Temperature = 18.60 


SPATIAL FREQUENCY* Boas Hate Cx 
(cycles/radian) (cm) (cm) (%) 
a0 Oe Oo ..05 URAC ace (0 G16) 7 0.20 2+ 20307 75.0 +e 
IS ole tel ce 2 1.40 OeZ7Z 72.8 foe 
2540 £045 1.40 0.2.6 68.7 Fez 
bi ciery 0 Mae came a>, 1.40 OrZe 66.7 ae 
a6 O bat e283 1.40 Orse 57.3 £2 
B64. Ome wlaae 1.40 0.44 52.2 fee 
Sg Osan Po 1.40 0.46 50.5 + 92eee 
7c 8 ace Any, o: 1.40 0.48 48.9 Lez 
6-35.07 22 3..6 1.40 OS. 0) 47.4 £ Zee 
LOO« 25 1.40 Osis: 41.4 £9235 
IW E38 1.40 0760 40.0, +15 
ZS <t 24 140 ONGZ 38.6 ee 
b3i6 cto es I. 40 0.64 37.3 toe 
LA Ft.55 eo) 0 21516 35.9 fee 
LG. EtG 1.40 0270 33.3 2 
TSS. £25 1.40 0.74 30.8 +e 
Z00 +9 1.40 0.84 25.0 £3 
Za So Ge 7 1.40 0.86 23.9 + eae 
Z22. sees 1.40 0.668 22.8 £0re 
250 + 9 1.40 0.94 19.7 £ Ome 
241). ae LO 1.40 0.96 18.6 2+ 30Ree 
256-2 9 le 40 0.98 L7.6 2Gee 
300) 26; 123 1.40 LO 12.0 £30Rs 
ile), “ashe 1.40 Vr 19 7.69 + ORs 
S50 7 aoe ae Sie) Lis 7.09 +30 
S506 i -5 Ies@ NS Le 6.56 £ Ue 
375 ot 216 O65 0-610 4.00 + OR 
400 + 17 0.60 0.60 O09 


x @f = (1/2w)dR + (R/2w2)dw 
Where df is estimated frequency error, w is bar width, R is range 
target, dw = 0.0002 m, and GR = 0. 0lem: 


** dC = [2(Hpay-Hmin)/(HmaxtHmin)*] GH + (df/f)C 
Where dC is Be oimewed radiation Contrase -crre.. 


SZ 


ieee 7 


NORMALIZED RADIATION CONTRAST AND MTF 


SPATIAL FREQUENCY NORMALIZED C MTF* 
(cycles/radian) (%) (%) 
200i tt 0.05 hO@at 3 99.7 + 
ieee Os 2 Sei t  2i7 Oi ot 
2.5 aid Ore Seon st 2.409 ‘eite) 4 10) eas 
Bore +t 0.9 Semone 3. 3 Bin. 3. -£ 
46.0 + 1.3 Toast 3.1 6922-1 
Bie oe 2 | 1.68 Semen: 2.7 es ON ae 
SeseeOete 2.2 Ove dee. 3 1 eae J Se 
eo 2.6 Gonz 3. 2 5S. 402 
Son) ot 3:,6 Gira: 935.5 Bice Set 
POO at -5 Sowzet 36.3 TA Ages a 
A este 5 Sioa 2. O 44.3 £ 
e254 Boe. 241 aL Js) ac 
ieee 3 ASeweet 1.6 390m 
iaset 5 ay emope ts. 2.1: STG 
Loi + 6 G42 -2 2.0 34.973 
Pest 5 Aneta S 3 2eaess ee 
Z200 +.9 Soom la. 2.02. 
ales) st. 7 Sloe Ll. 3 Ds, se 
Ze2 t 5 SOm4e 2 0% 9 2 ae) Ge 
Z50) + 9 Zowe t 1.2 206 
27 £ 10 Zoe ee lL 19-6. + 
286 + 9 Zoot 0.9 18.5 + 
S00. + 13 Mee. t+ O.. 8 26. 
Seat, 1) Ores tO. 4 Grae get 
3356+ 12 O245 7+ 0.40 (G2 se 
S505+ 15 Sato t 0.40 Gero at 
So) = 16 Seo t+ -O)..30 QO oe 
400 + 17 O20: 


x AMTF = dc 
Where AMTF is the estimated MTF error. 
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SPATIAL SC REOGENGY 
(cycles/radian) 


* 


if 
ON 
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GaMRTD 
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= (4.00 X 107-4) [d£/MTF + (£/MTE) 
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Where AMRTD is the estimated MRTD error. 


WNYRRrF WON MOP WW WNHRRRRPPRP OUP BP WNF FR OF 
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NO 
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dMTF ] 
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TABLE 9 
Revlon TEMPERATURE DIFFERENCE (°C) 
(THERMOVISION TEMPERATURES - THERMISTOR TEMPERATURES ) 


DATE 6 May 9 May 11 May 13 May 
MAERMISTOR 2 - 0.4 eeOiae ap ANG, ~ (as 
MAE RMISTOR Z - 1.0 = Oe dt sO 9 — 10.6% 
AP RMT STOR 3 Se, =, 2, 0 68 = 105 fe 
THERMISTOR 4 S057 - 0.4 + 0.7 - 1.4 
MAERMISTOR 5 - 2.9 = NO) te sp Ol SS, a tae 
THERMISTOR 6 SSCs = (52 =) 0h 8 =O 5 
THERMISTOR 7 = 6.3 * * * 

THERMISTOR 8 Se Or * * * 


* Not available. 
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